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FIRESTORM  ANALYSIS 


1.0  INTRODUCTION 

1. 1  Objectives 

The  following  general  objectives  were  specified  for  this  study: 

a)  Compute  the  radial  inrush  wind  (vw)  asso¬ 
ciated  with  a  mass  fire  as  a  function  of  time 
and  various  fuel  bed  parameters, 

b)  Investigate  the  sensitivity  of  vw  to  varia¬ 
tions  in  certain  of  the  fuel  bed  parameters, 
and 

c)  Conduct  an  analytical  study  of  the  case  in 
which  fire  spread  is  accelerated  by  wind 
and/or  flame  interaction  until  coalescence. 

1. 2  Background 

The  basic  point  of  departure  for  this  research  was  Ref.  1,  in  which 
a  hypothesis  was  developed  relating  radial  inrush  wind  velocity  at  the  edge 
of  a  potential  firestorm  area  to  the  energy  release  rate  of  fires  within  the 
area  and  to  the  size  of  the  area  itself.  The  value  of  the  constant  in  the  equa¬ 
tion  relating  these  parameters  was  determined  by  "fitting"  an  equation 
derived  from  physical  considerations  to  estimated  "worst  case"  conditions 
during  the  Hamburg  firestorm.  One  degree  of  freedom  in  the  "constant" 
term  allowed  this  fit  to  be  made. 

In  Appendix  IB,  the  equation  of  Ref.  1  is  expanded  to  express  the 
energy  release  rate  of  the  fires  in  terms  of  the  following  parameters: 


n 


the  number  of  initial  fires, 


FL  the  fuel  loading, 

E.  j  the  energy  content  of  the  fuel, 

A  the  area  of  a  single  fire  which  is  in  the 

y  x 

violent  burning  stage  at  a  given  time, 

the  area  of  a  single  fire  which  is  in  the 
residual  burning  stage*  at  a  given  time, 

E.  .  ,  the  fraction  of  Efllc,i  released  during  the 
trac,  viol  .  ,  iuei *  0 

violent  burning  stage  , 

E  the  fraction  of  Efuej  released  during  the 

rac,  resi  residual  burning  stage*. 

At  the  duration  of  the  violent  burning  period,* 

and 

AtD  the  duration  of  the  residual  burning  period*, 

ft 


In  Appendices  IA  and  IC,  equations  are  derived  which,  for  certain 

assumptions,  represent  A  and  A  as  a  function  of  time,  spread  rate, 

SY  o  ft 

and  other  parameters  for  three  different  spread  conditions. 


1.  3  Method  of  Calculation 

A  computer  program  to  calculate  v  was  written  based  upon  Eq. 
(IA-1)  and  Tables  IA-I  and  IA -H .  The  input  values  used  are  listed  in  Sec¬ 
tion  IA.  4.  The  results  of  these  calculations  are  presented  graphically  in 


The  definitions  of  "violent"  and  "residual"  burning  regimes  are  based 
upon  Ref.  2. 


•  2- 


Appendix  IV.  A  second  part  of  the  program  calculates  the  partial  derivative 
of  v^  with  respect  to  seven  selected  parameters  based  upon  the  rationale 
of  Appendix  II  and  using  the  inputs  listed  in  Table  II-I.  The  results  of  this 
series  of  calculations  is  presented  in  Appendix  V.  A  listing  of  this  pro¬ 
gram,  its  input  requirements,  and  a  sample  output  page  from  each  part  are 
contained  in  Appendix  III . 

The  accelerated  spread  study  considers  that  radial  spread  is  accel¬ 
erated  by  interaction  effects  (after  adjacent  fires  reach  a  certain  minimum 
separation  distance).  Four  initial  spread  modes  are  investigated:  a  con¬ 
stant  rate  of  radial  increase,  a  constant  rate  of  area  increase,  and  two  log- 
arithmic  rates  of  radial  spread. 

The  computer  program  iisted  in  Appendix  VII  calculates  spread  rate, 
fire  radius,  and  distance  to  the  point  at  which  flame  fronts  merge  as  a 
function  of  time  for  each  initial  spread  mode.  This  program  is  based  upon 
the  equations  of  Appendix  VI .  Appendix  VIII  contains  eight  graphs  showing 
radial  spread  rate  as  a  function  of  time  and  of  fire  radius  for  each  of  the 
initial  spread  modes,  with  the  point  at  which  acceleration  begins  indicated 
on  each  curve. 

1. 4  Limitations 

Inherent  in  the  rationales  and  the  results  of  thi  study  are,  of  course, 
the  approximations  and  assumptions  made  in  deriving  the  heuristic  model 

* 
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from  Ref.  1  and  the  modifications  to  it  described  in  Section  1.2.  In  partic¬ 
ular,  the  "worst  case"  conditions  in  Hamburg,  from  which  the  value  of  the 
constant  was  obtained,  corresponds  to  100%  simultaneous  involvement  of 
all  structures  in  an  urban  area  during  the  violent  burning  period  of  the  ini¬ 
tial  fires.  This  study  considers  a  continuous  wildland  fuel  bed  with  some¬ 
what  different  burning  characteristics.  The  degree  of  correspondence 
between  these  two  cases  is  not,  as  yet,  firmly  established. 

All  of  the  v  computations  are  made  using  average  values  for  the 
rate  of  energy  release  term,  as  reflected  by  the  parameters  listed  on  pages 
IA-1  and  1A-2. 

Also,  the  type  and  rate  of  fire  spread  which  occur  in  a  fuel  bed  of 
the  sort  considered  here  will  probably  r  quire  considerable  further  inves¬ 
tigation  before  the  functional  dependencies  and  appropriate  parameter  values 
are  clearly  defined.  The  preliminary  work  on  accelerated  spread  in  this 
study  is  one  course  of  investigation  that  may  yield  a  better  understanding 
of  these  relationships. 

The  spread  rate  relations  assumed  are  approximations,  and  in  the 
real  world  probably  represent  rough  average  values  as  they  develop  as  a 
function  of  time. 

This  entire  study  may  be  viewed  as  a  project  that  has  generated  a 
number  of  rough  hypotheses  which  may  now  be  tested  and  assessed  as  to 
validity  and  utility. 
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2. 0  RESULTS 


2. 1  Radial  Inrush  Wind  Velocity  Curves 

The  45  families  of  curves  in  Appendix  IV  represent  as  a  func¬ 
tion  of  n  for  ten  values  of  time  after  ignition.  A  family  of  curves  was  con¬ 
structed  for  each  combination  of  values  of  3  spread  cases,  3  fuel  loadings, 
and  5  area  sizes.  Each  family  is  bounded  above  by  the  coalescence  locus, 
as  discussed  in  Appendix  ID. 

In  general,  an  increase  in  the  number  of  initial  ignitions  decreases 

the  time  at  which  the  maximum  value  of  v  (v^  )  occurs,  as  would  be 

\v  w 

expected.  For  Case  1  (constant  area  increase),  v^  (as  described  by  the 
coalescence  locus)  increases  considerably  as  a  function  of  n.  On  the  other 
hand,  for  Case  2  (constant  rate  of  perimeter  or  radial  increase),  the  coa¬ 
lescence  locus  is  much  "flatter";  i.  e. ,  v^  is  fairly  constant  with  respect 
to  n.  This  simply  means  that,  as  n  increases,  the  time  to  coalescence  is 
less  for  spread  Case  1  than  for  spread  Case  2,  all  other  factors  being 
equivalent.  Also,  for  the  values  of  the  spread  constants  considered  in  this 
analysis,  v^  appears  to  occur  at  a  much  earlier  time  for  Case  1  than  for 
Case  2,  which  is  intuitively  expected. 


2.  2  Sensitivity  Analysis 


Appendix  V  contains  seven  families  of  curves  which  represent 


at  ten  times  after  ignition  as  a  function  of  jj.  with  the  other  parameters 
pr  (j  i  i)  held  at  their  "representative  values.  "  These  curves  are  drawn 


from  the  output  of  Part  II  of  the  computer  program  listed  in  Appendix  IE. 
This  part  of  the  program  and  its  inputs  are  based  upon  the  rationale  des¬ 
cribed  in  Appendix  II. 

9v 

w 

The  partial  der  ivative  — —  indicates  the  rate  of  change  of  v 

9u .  w 

9v 
w 

with  respect  to  ju.  .  The  sign  of  — —  indicates  the  direction  of  change, 

9v  1  d^i 

w 

e.  g.  ,  if  *r —  >0,  an  increase  in  the  value  of  y.  will  result  in  an  increase 

dyi  9v  1 

w 

in  the  value  of  v  .  Conversely,  if  — —  <  0,  then  an  increase  in  n-  wil 
w  J  9 ju,  .  «  . 


yield  a  decrease  in  v^ .  The  magnitude 


9v 


w 


9  Hi 


serves  as  an  index  of  the 


rate  of  change  for  various  points  on  the  family  of  curves  associated  with  a 

9v 

w 

particular  Table  2.2-1  lists  the  sign  of  — —  and  the  behavior  of 

l  au . 

l 


9v 
_ w 

9wi 


for  the  seven  parameters  considered. 


2.  3  Accelerated  Spread  Study 

Appendix  VI  considers  the  case  in  which  radial  spread  rate  is  accel¬ 
erated  by  interaction  effects  when  adjacent  fires  are  sufficiently  close. 

This  case  is  investigated  for  four  initial  spread  modes;  two  logarithmic  rates 
% 

of  radial  increase,  the  constant  rate  of  radial  spread  (M  =  0.  01  ft/sec) 

used  for  Case  2b  of  the  v  calculations,  and  the  constant  rate  of  area 

w 

2 

spread  (K,  =  5.  0  ft  / sec)  used  for  Case  1  of  the  v  calculations. 

1  w 
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from  the  output  of  Part  II  of  the  computer  program  listed  in  Appendix  III. 


This  part  of  the  program  and  its  inputs  are  based  upon  the  rationale  des¬ 
cribed  in  Appendix  II. 

av 

W 

The  partial  derivative  — —  indicates  the  rate  of  change  of  v 

8“i  8v 
w 

with  respect  to  /u.  .  The  sign  of  — —  indicates  the  direction  of  change, 

3v  1  0Mi 

w 

e.  g. ,  if  -  >0,  an  increase  in  the  value  of  y.  will  result  in  an  increase 

9^i  8v  1 

w 

in  the  value  of  v  .  Conversely,  if  — —  <  0,  then  an  increase  in  u .  will 
w  J  dn.  .  _  i 


yield  a  decrease  in  v  .  The  magnitude 


9v 


w 


serves  as  an  index  of  the 


rate  of  change  for  various  points  on  the  family  of  curves  associated  with  a 


3v 


w 


particular  Table  2.2-1  lists  the  sign  of  ^ 


and  the  behavior  of 


8v 


w 


for  the  seven  parameters  considered. 


2.  3  Accelerated  Spread  Study 

Appendix  VI  considers  the  case  in  which  radial  spread  rate  is  accel¬ 
erated  by  interaction  effects  when  adjacent  fires  are  sufficiently  close. 

This  case  is  investigated  for  four  initial  spread  modes;  two  logarithmic  rates 

a!c 

of  radial  increase,  the  constant  rate  of  radial  spread  (M  =  0.  01  ft/sec) 

used  for  Case  2b  of  the  v  calculations,  and  the  constant  rate  of  area 

w 

2 

spread  (K,  =  5.  0  ft  / sec)  used  for  Case  1  of  the  v  calculations. 
i  w 
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Equations  are  developed  for  radial  spread  rate,  fire  radius,  and 
distance  to  merging  as  a  function  of  time  for  each  initial  spread  mode.  The 
computer  program  listed  in  Appendix  VII  is  based  upon  these  equations. 

The  eight  graphs  in  Appendix  VIII  represent  radial  spread  rate  as  a 
function  of  time  and  of  fire  radius  for  each  initial  spread  mode. 
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APPENDIX  IA 


r 

i  * 


RATIONALE  FOR  COMPUTATION  OF 
RADIAL  INRUSH  WIND  VELOCITY 


i 

I 

I 

I 

r 


IA.  1  Basic  Equation  and  Parameter  Definitions 

* 

The  basic  equation  is 


1  88 

50  ] 

1/3 

[(nKFLKEfuei)] 

fA  (r,t,At  )*E  . 

SV  V  frac,  viol 

V  ~ 

w 

[1000 

1089 

Al/8 

ASR(r,t’  AV  AV 


-.1/3 


Trac,  resid 


At 


R 


(IA.  1) 


Where 

v  is  the  radial  inrush  wind  velocity  at  the  periphery 

W  of  a  circular  area  A  within  which  fires  are  burn¬ 

ing,  in  miles  per  hour 

n  is  the  number  of  initial  fires.  [  n^^-  (A)  ,  where 

n  ,  is  the  number  of  initial  fires  per  square  milel 
d 

FL  is  the  fuel  loading,  in  tons  per  acre 

E„  ,  is  the  caloric  content  of  the  fuel  in  Btu/lb 

fuel 

A  is  the  total  area  within  which  fires  are  started,  in 

square  miles 

Agy  is  the  area  which  is  in  the  violent  burning  regime 

at  time  t,  in  square  feet  ■ 

i 


This  equation  is  based  on  work  done  in  Ref.  1.  A  derivation  of  the  form 
used  in  this  study  from  those  in  Ref.  1  is  given  in  Appendix  IB. 


IA-1 


r 


is  the  spread  rate;  if  Ag  denotes  the  total  area, 
in  square  feet,  over  which  the  fire  has  spread 
at  time  t , 


r 


dt 


t  is  the  time  after  ignition,  in  seconds 

Aty.  is  the  duration  of  the  violent  burning  period,  in 

seconds 


Efrac  vj0}  is  the  fraction  of  the  total  available  energy  of  the 
’  fuel  which  is  released  during  the  violent  burning 

period 


is  the  area  which  is  in  the  residual  burning 
regime  at  time  t,  in  square  feet 

At^  is  the  dura+ion  of  tne  residual  burning  period,  in 

seconds 


E„  .  ,  is  the  fraction  of  the  total  available  energy  of  the 

i r*3.c  resid 

'  ‘  ’’ .  fuel  which  is  released  during  the  residual  burn¬ 

ing  regime 


IA.  2  Assumptions 


A.  All  fires  have  negligible  initial  areas;  i.  e.  ,  they  are  effec¬ 
tively  "point  fires". 

B.  Fire  spread  starts  immediately  after  initiation  of  any  given 
"point  fire",  and  each  fire  spreads  in  the  same  manner  at 
the  same  rate  symmetrically  around  the  "point  fire"  until 
coalescence. 

C.  All  fires  are  started  simultaneously  from  a  uniform  distri¬ 
bution  of  ignitions  and  burn  in  identical  fuel  beds. 


(IA.  2) 


IA-2 


D.  Three  rates  of  area  spread  with  time  will  be  considered: 


dA0 

_ _ O  _ 

dt 

dA 

ri  ■  ki  • 

which  represents  a  constant 
rate  of  area  growth, 

o  _ 

dt 

r2  =  Bt' 

which  represents  a  constant 
rate  of  perimeter  increase. 

dA 

_ O 

dt 

r3  =  Ct, 

which  represents  a  constant 
rate  of  radial  spread. 

The  expressions  for  B  and  C  are  derived  in  Appendix  IC. 


IA.  3  Procedure 

At  any  time  t ,  the  area  into  which  a  given  fire  has  spread,  Ag ,  is 
composed,  ideally,  of  three  regions: 


1) 

2) 

3) 


The  violently  burning  part,  Agy  , 

The  residually  burning  part,  A  ,  and 

on 

The  burned- out  portion,  AeTJ^  ,  so  that 

oBU 


AS(t)  “  ASV<t}  +  ASR(t)  +  ASBO(t) 


From  the  definitions,  it  is  apparent  that 


A  (t)  =  A  (t)  -  A  (t  -  Aty) 

ASR^^  "  AS(t'AV  ’  As[t_  ^Atv+ AtR^ 
ASBO(t)  =  As[ t_  (Aty+  AtR^ 


(I A.  3) 


(1A.  4) 
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Taking  r(t)  as  the  spread  rate,  one  may  write  the  following  rela¬ 


tions: 


SBO 


0  if  t  <  Aty  +  AtR 


(Aty  +  AtR)]dt  for  t  >  Aty  +  AtR 


(IA.  5) 


AtV  +  AtR 


SR 


0  if  t  <  Aty 


r 


(IA.  6) 


r(t  -  Aty)  dt  -  AgBO  for  t  >  Aty 


At, 


and 


SV 


r 


r(t)dt  -  Asr  -  Asbo 


(IA.  7) 


These  expressions  have  been  evaluated  and  the  results  are  given  in 

Tables  IA-1  and  IA-II,  which  present  formulae  for  A  ,  ,  A  ,  and  A  as 

oV  bn  u  dU 

a  function  of  time  for  the  two  basic  spread  modes  given  in  paragraph  IA.  2D 
above. 


Case  1. 


Constant  area  increase. 


Table  IA  - 1 


0  <  t  <  Aty 

Aty  <  t  <  Atv  +  AtR 

1  >  AtV  +  AtR 

Asv 

rt 

r  •  Aty 

r  •  Aty 

asr 

-0- 

r  •  (t  -  Aty) 

r-AtR 

asbo 

-0- 

-0- 

r  *  [t  -  (Aty  + AtR>] 

In  this  case 


nr  =  ~  X  (5280)2 


(IA.  9) 


where  is  the  number  of  initial  fires  necessary  to  achieve  coalescence 
in  area  A  by  time  t . 


dA 

Case  2.  r  =  — r~  =  K„t 

- —  dt  i 


(I  A.  10) 


2 

(2a)  Kg  =  B  (constant  rate  of  perimeter  increase)  B=L  I  (2*) 
where  L  is  rate  of  change  of  perimeter,  i;i  feet 
per  second. 


2 

(2b)  Kg  “  C  (constant  rate  of  radial  increase)  C  =  2irM 

where  M  is  rate  of  change  of  radius,  in  feet  per 
second. 


The  derivation  of  the  coalescer.ce  relations  is  given  in  Appendix  ID. 
The  derivation  of  these  relationships  :  s  given  in  Appendix  IC. 
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4200  sec 

2, 

5.0  ft  /sec 
.01  ft/sec 
.1  ft/sec 
0.4 
0.6 

2,  4,  6,  8,  16,  24,  32,  48,  64,  80  minutes 

Note  that  the  total  number  of  fires  burning  in  a  given  area  is  repre¬ 
sented  by  n  =  n^  •  A  . 

IA.  5  Method  of  Calculation 

A  computer  program  was  written  based  upon  the  above  equations. 

Using  the  inputs  listed  above,  it  calculates  v  as  a  function  of  time  and  n 

w 

for  each  combination  of  the  other  parameters.  A  listing  of  this  program, 
its  input  format  requirements,  and  a  sample  output  page  is  included  in 
Appendix  III.  The  output  is  grouped  so  that: 

1)  All  pages  pertaining  to  a  given  spread  rate 
are  together. 

2)  Within  a  given  case,  "’ll  pages  with  a  specific 
FL  are  together. 

3)  Each  page  of  the  output  has  all  the  points  for 
a  specific  area  size,  given  the  case  and  FL. 

The  results  of  these  calculations  are  presented  graphically  in 

Appendix  IV. 


At 


R 


K 


1 

M 


E 

frac, viol 

^frac,  resid 
t 
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APPENDIX  IB 


Equation  (29)  of  Ref.  1  gives 


(88  X  10"3) 


n 


l 


dE, 


1 1/3 


T(i) 


dt 


L  i=l 


w 


(IB.  1) 


where  is  the  radial  inrush  wind  velocity  at  the  periphery  of  area  A,  in 

miles  per  hour,  n  is  the  number  of  initial  fires  started  in  an  area  A,  and 
dET(i)  , 


dt 


is  the  energy  release  rate  for  a  single  fire,  in  Btu/sec  at  time  t. 


Assuming  identical  simultaneously  ignited  fires  in  a  continuous  fuel  bed  with 
uniform  characteristics,  the  total  amount  of  energy  released  during  the  vio¬ 
lent  burning  period  is  given  by  the  product  of:  - 


DERIVATION  OF  THE  BASIC  EQUATION 


1)  the  total  number  of  fires, 

2)  the  fuel  loading  (lbs /ft3), 

3)  the  caloric  content  of  the  available  fuel 
(Btu/lb), 

2 

4)  the  area  (ft  )  which  is  burning  violently  at 
the  time  in  question,  and 

5)  the  fraction  of  the  total  energy  of  the  fuel  which 
is  released  during  the  violent  burning  period. 

The  average  energy  release  rate  during  the  violent  burning  regime  is  the 
total  energy  released  during  that  time  divided  by  the  length  of  the  violent 
burning  period,  in  seconds.  Symbolically, 


IB-1 


I 


dE, 


T(i) 


dt 


>-  i=  1  J  viol 


(n)(FL)(E  )(A  )(E  .  ) 

_ fuel  5V  frac,  viol 

AtTr 


(IB.  2) 


Similarly,  the  energy  release  rate  during  the  residual  burning  period  may 
be  represented  by 


r  n 


I 

L  i=l 


dE, 


T(i) 


dt 


M<FL><Efuel)(ASR)(Efrac.  resid1 
At„ 


(IB.  3) 


resid 


If  fuel  loading  (FL)  is  expressed  in  tons  per  acre,  the  expressions 
on  the  right  in  Eqs.  (IB.  2)  and  (IB.  3)  must  be  multiplied  by  the  conversion 
factor  (50/1089). 

The  total  energy  release  rate  is  the  sum  of  the  release  rate  for  those 
portions  of  the  fuel  bed  in  the  violent  burning  regime  and  for  the  portions  in 
the  residual  regime,  or 


V  dET(i) 

[r  dET(i)l 

r  n  -i 

V  d  T(i) 

Z  dt 

Z  dt 

T 

Z  dt 

Li=i  J 

total 

Li=l  J 

viol 

Li=i  J 

resid 

Substituting  Eqs.  (IB.  2),  (IB.  3),  and  (IB.  4)  into  Eq.  (IB.  1),  using  the  con¬ 
version  factor,  and  simplifying  yields 


88  | 

|_50J 

l>/3  |(n)(FL)(Efuel)  | 

1 1/3 

^SV  Efrac,viol  ,  ^SR  Efrac,  resid 

ilOOO/ 

1 1089] 

1  he 

L  A*V  **R  J 

which  is  the  desired  result.  In  this  expression  is  a  complex  function 

of  r,  t,  and  At^. ,  and  Ag^  is  a  complex  function  of  r,  t,  At^,  and  At^ 


1/3 

(IB.  5) 


IB- 2 


r^i 


APPENDIX  IC 


DEVELOPMENT  OF  SPREAD  RATE  RELATIONSHIPS 

The  expressions  for  B  and  C ,  as  defined  under  assumption  (D), 
Section  IA.  2,  may  be  derived  as  follows: 


a)  Let  the  perimeter  of  a  fire  be  approximated  by  P  =  2?rrc  , 

where  r  is  the  radius  of  the  circular  fire  front  and  is  a  function  of  time, 
c 

If  the  rate  of  perimeter  growth  is  constant,  then 

dP 


dt 


=  L 


(IC.  I) 


I 


r 


!  . 


where  L  is  the  constant  rate  of  perimeter  increase  in  units  of  length  per 
unit  of  time. 

From  Eq.  (IC.  1), 


and 


Now, 


P  =  Lt  =  2nr 


Lt 

2tt 


.  2  LV 

AS  =  "c  =  — 


(IC.  2) 


(IC.  3) 


(IC.  4) 


Differentiation  of  Eq.  (IC.  4)  yields 


dAc 

_ k. 

dt 


L2t 

2  it 


r„  =  Bt 


(IC.  5) 


IC-1 


,v.y .  if 


Solving  for  B  yields 


B  = 


L 

27T 


(ic,  6) 


b)  Similarly,  if  r^  again  represents  the  radius  of  a  circular  fire 
front  and  M  represents  a  constant  rate  of  radial  spread  in  units  of  length 
per  unit  of  time,  one  may  write 


1. 

r 

i. 

[ 

i 

[ 


so  that 


and 


Hence 


dr 

_ c 

dt 


=  M 


r  =  Mt 
c 


2  2  2 
A„  =  ?rr  =  ttM  t 
S  c 


dA 


=  27rM2t  =  r  =  Ct 

dx  u 


and  therefore 


C  =  2ttM 


(IC.  7) 


(IC.  8) 


(IC.  9) 


(IC.  10)  | 


(IC.  11)  | 
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DERIVATION  OP  COALESCENCE  EQUATIONS 


The  v  versus  n  curves  are  bounded  by  the  "coalescence  locus", 
w 

which  may  be  defined  as  the  number  of  initial  fires  necessary  to  achieve 
* 

coalescence  by  a  specified  time. 

Assuming  a  uniform  density  of  (n)  ignitions,  the  area  in  square  feet 
available  to  a  single  fire  for  spread  before  coalescence  is 


a  _  A  X  (5280)2 
AS  n 

Also, 


Ag(t)  =  ^  r(t)  dt 


(ID.  1) 


(ID.  2) 


If  r(t)  =  Kx  (Case  1), 

Ag(t)  =  Kjt 


(ID.  3) 


Equating  (ID.  1)  and  (ID.  3)  yields,  for  Case  1, 


n 

c 


A  X  (5280) 
Ki‘ 


(ID.  4) 


At  coalescence  all  of  the  individual  fires  have  merged  into  a  single  large 
fire,  and  it  is  assumed  that  no  regions  within  the  area  remain  unburned. 


ID-1 


where  nc  represents  the  approximate  number  of  initial  fires  required  for 
coalescence  by  time  t.  The  approximation  results  from  geometric  assump¬ 
tions  inherent  in  Eq.  (ID.  4). 

In  Case  2, 


r  =  K2t 


and  Eq.  (ID.  2)  becomes 


As(.) 


■  I‘ 

w  n 


V dt  =  1  S'2 


(ID.  5) 


(ID.  6) 


Equations  (ID.  1)  and  (ID.  6)  yield  the  approximation 


n 

c 


2A  X  (5280) 


2 


K2‘ 


.2 


for  spread  Case  2. 


(ID.  7) 
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APPENDIX  II 

RATIONALE  FOR  SENSITIVITY  ANALYSIS 


II.  1  Purpose 

To  examine  the  sensitivity  of  the  radial  inrush  wind  velocity,  as 


calculated  in  Appendix  IA,  to  variations  in  seven  selected  parameters, 

0  V 
w 

expressions  for  the  partial  derivatives  — —  have  been  developed,  where 


A 


u3=FU 


u=  E. 

4  fuel 


U  =e 

5  frac,viol’ 


^6  Efrac,resid’ 


and  p7=  r  . 


II.  2  Partial  Derivatives 

Formulae  for  the  partial  derivatives  with  respect  to  selected 
parameters  have  been  derived  as  follows: 

For  convenience,  let  us  represent  Eq.  (IA.  1)  as 


vw  ’  V  Q2  ’  Q3 


where 


Q, 


88 

50  1 

1000 

1  089  J 

1/3 


Q_ 


[  (n)(FL)(Efuei)  ] 

a  1/6 


1/3 


and 


Q,  =  { 


ASv(r,t,Aty)’  Efrac,  viol  +  ASR(r AV  AV  ‘  Efrac,  resid 


AtT 


At 


1/3 


R 


II- 1 


(II.  1) 


(II.  2) 

(II.  3) 


(II.  4) 


V 


/ 


Then 


[•KLXE  '] 

Q1Q3  ,  2/3'V^S 
3n  A 


(II.  5) 


Q1Q3[(n)(FL)(Efuel)]' 
*  ,7/6 


-6  A 


(II.  6) 


3(FL) 


SQ3[(n)(Efuel^ 
3(FL)2/3  A1/6 


(II.  7) 


8(IW 


Q1Q3[(n)(FL)]1/‘ 

3(E.  /V* 

fuel 


(II.  8) 


a(E,  .  ,) 
frac,  viol 


Q1Q2  lAS  V(r't,AtV)  ASR^r,t*  AtV'  AtR* 

3(Q3)2  AtV  AtR 


(II.  9) 


^Efrac,  resid^ 


Q1Q2  /  ASR(r,t’AtV'AtR)  Asv(r.t.Aty) 

3<Q3)2  i  AtR  AtV 


(II.  10) 


Note  that 


9(Ef  .  J 
frac,  viol 


9(Efrac,  resid) 


(II.  ID 


For  Case  1 ,  if  t  <  Aty  , 


Q1Q2  1  ‘  Efrac,  viol 
3(Q3)2  AtV 


(II.  12) 


If  Aty  <  t  <_  Aty  +  AtR  , 


dv 

_ w 

9r 


SQ2 

3<Q3)2  L 


(t  Aty)  ^frac  resjd 

pp  -f.  .  ■  ..  -  — ■ . A  ■  —  —  ■ 

frac, viol  At„ 


and  if  t  >  At,T  +  At 


Aty  +  AtR  , 


3v 

_ w 

9r 


Q1Q2 

3(Q3)2 


E  E 

\  frac,  viol  frac,  resid 


Since 


^frac, viol  ^frac,  resid 


then 


9v 

w 

3r 


SQ2 

3(Q3)2 


For  Case  2,  if  t  <  Aty  , 


8v 

_ w 

8r 


S  Q2 

6<q3)2 


t  •  E 


frac,  viol 


At 


V 


If  Aty  <  t  <  Atv  +  AtH  , 


a  v 

w 

ar 


Q1  Q2 


2  -  (At 


6(Q3)“  L 


./t)  E'-n  v 

V  I  fiiic.v 


iol 


t  -  2Aty  +  (Aty It) 


At 


H 


!•: 


frac,  resid 


(II.  13) 


(II.  14) 


(II.  1 5) 


(n.  16) 


(II.  17) 


.(II.  18) 


II- 3 


and  if 


i  >  Afy  +  AtR  , 


av 

_ w 

ar 


Q1  Q2 

6(Q3)2 


r 

2(At  )  +  At  " 

<  2  -  (Aty/t) 

E  + 

frac,  viol 

n  v  n 

t 

E  .  ,  >  . 

frac,  resio 


II.  3  Input  Parameters 

Calculations  have  been  made  of 


av 


w 


i  =  l  , 


7 ,  using 


fixed  representative  values  of  (j  4  i)  and  varying  p.  across  ;  range. 
Representative  values,  ranges,  and  points  of  calculation  are  listed  in  Table 


II-I. 


II.  4  Method  of  Calculations 

A  computer  program  was  written  based  upon  Eqs.  (II.  2),  •  •  •  ,  (II.  19) 

and  run  using  the  input  parameters  as  defined  in  Table  II-I.  A  listing  of 

this  program  and  a  sample  of  its  output  is  included  in  Appendix  III. 

Appendix  V  contains  seven  sets  of  sensitivity  curves  based  upon  the 

9v 

W 

output  of  this  program.  Each  set  represents  plotted  as  a  function  ot 

i 

p.  for  each  time  value. 


(II.  19) 
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TABLE  l!  -  1 


INPUT  VALUES  KOR  SENSITIVITY  ANALYSIS 


Representative 

Parameter 

Value 

Range 

Points 

nd 

4500 

2K-  10K 

2K,  3K,  4K,  4.5K, 
5K,  6K,  7K,  and 
10K 

FL 

15.0  T/Ac 

5-30  T/Ac 

5,  15.0,  30  T/Ac 

E.  , 
fuel 

7000  Btu/lb 

6000  -  9000 

6K,  7K,  8K,  9K 

A 

15.0  mi2 

2 

1 .0  -  25.0  mi 

1.0,  5.0,  10.0 

15.0,  25.0  mi2 

t 

-- 

(Same  as 

in  Section  IA,  4) 

> 

<1 

360  sec 

-- 

r 

(Case  2-b) 

M 

.01  ft/sec 

.00  5 -.2  ft /sec 

.005,  .01,  .015,  .2 

L 

.1  ft/sec 

-- 

-- 

Ki 

5.0  ft-/see 

Asv 

(Use  Tabic  IA-II) 

"Trac,  viol 

0.4 

0.2  -  0.8 

0.2,  0.4,  0.6,  0.8 

aik 

4200  sec 

-- 

-- 

ASli 

(Use  Table  1A -11) 

Trac,  res  id 

0.6 

0.2  -  0.8 

0.2,  0.4,  0.6,  0.8 
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APPENDIX  III 


COMPUTER  PROGRAM  FOR  INRUSH  WIND  AND 
SENSITIVITY  CALCULATIONS 


III.  1  General 

The  following  program  is  in  two  parts.  Part  I  calculates  v  as  a 

function  of  time  and  n  for  each  combination  of  the  other  input  variables  as 

9  v 

w 

defined  in  Section  IA.  4  .  Part  II  calculates  — —  (i  =  1  ,  ■  •  •  ,  7  )  as  a 
function  of  and  time,  with  p_.  (j  4  i)  held  constant,  as  discussed  in  Sec¬ 
tion  II.  4.  The  inputs  for  Part  II  are  listed  in  Table  II-I.  The  program  is 
written  in  FORTRAN  IV  for  an  IBM  System/ 360  Model  40G  digital  computer. 


HI.  2  Data  Cards 

Thirteen  data  cards  are  required.  The  first  six  data  cards  contain 
the  input  data  for  Part  I  of  the  program.  Card  1  has  format  8F10.  4  and 
contains  the  following  variables: 


Columns 

Parameter 

Program 

Symbol 

Units 

1-10 

Lfucl 

i ;  i 

Btu/lb 

11-20 

Atv 

OTV 

seconds 

21-30 

AtR 

PTR 

seconds 

31-40 

*'frac,  viol 

KFV 

number 

41-50 

*Jfrae,  resit! 

KFR 

number 

III- 1 


Card  2  has  format  3F10.  6  and  contains,  in  cols.  1-10,  the  rate  of 
area  spread  in  square  feet  per  second;  in  cols.  11-20,  the  rate  of  perimeter 
increase  in  feet  per  second;  and  in  cols.  21-30,  the  rate  of  radial  increase 
in  feet  per  second. 

Card  3  has  format  8F10.  4  and  contains  three  values  of  fuel  loading 
in  tons  per  acre. 

Card  4  has  format  8F10.  4  and  contains  five  values  of  area  in  square 

miles. 

Card  5  has  format  818  and  contains  eight  values  of  initial  ignition 
density  (n^)  in  number  per  square  mile. 

Cards  6  and  7  have  format  8F10.  4  and  contain  ten  values  of  time  in 
seconds. 

Card  8  has  format  818  and  is  a  flag  card.  If  one  wishes  to  perform 
the  sensitivity  calculations  of  Part  II,  an  integer  2  should  be  put  in  col.  8 . 

If  one  wishes  to  omit  the  sensitivity  calculations,  an  integer  1  should  be 
put  in  col.  8  . 

The  last  five  cards  contain  the  additional  inputs  needed  for  Part  II 
of  the  program.  Card  7  has  format  18,  6F10.  4  and  contains  the  represent¬ 
ative  values  of  the  seven  p.  arranged  in  the  following  manner: 
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Columns 

Parameter 

Program 

Symbol 

Units 

1-  8 

Hep.  value  of 

NDR 

number 

9-18 

Rep.  value  of  A 

AR 

square  miles 

19-28 

Rep.  value  of  FL 

FLR 

tons /acre 

29-38 

Rep.  value  of  E„  . 

fuel 

EFP 

Btu/lb 

39-48 

Rep.  value  of  E. 

frac,  viol 

EFVR 

number 

49-58 

Rep.  value  of  E,  .  , 

r  frac,  resid 

EFRR 

number 

59-68 

Rep.  value  of  M 

AMR 

feet/ second 

Cards  10-13  have  format  4F10.  5  and  contain  the  points  across  the 
range  of  the  four  p.  which  were  used  as  constant  in  Part  I. 


Program  No.  of 


Card 

Parameter 

Symbol 

Points 

Units 

8 

Efuel 

EF2(I) 

4 

Btu/lb 

9 

Efrac, viol 

EFV2(I) 

4 

number 

10 

Efrac, resid 

EFR2(I) 

4 

number 

11 

M  (radial  spread  rate) 

AM2(I) 

4 

feet/ second 

III.  3 

Program  Listing 
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DIMENSION  A(5) ,T( 10I,N(8) ,CN(8),FL( 3) ,NC( 10) ,ND(8) ,VW(10,8!,TM(10)  1 

DIMENSION  EF2v4), EFV? ( 4 ) , EFR 2( 4 ) , AM21 4 ) ,ASVR( 10) , ASRR ( 10 ) , 03 ( 10 )  2 

$i PVW ( R »  10  ) , NA ( 8 )  3 

C  PART  ONE  CALCULATES  THE  SET  OF  RADIAL  INDRAFT  WIND  VELOCITIES  4 

C  RESULTING  FROM  8  VALUES  OF  N,  THE  NUMBER  OF  INITIAL  FIRES»  AT  10  5 

C  GIVEN  TIMES  AFTER  IGNITION.  THIS  SET  IS  PRODUCED  FOR  COMBINATIONS  6 

C  OF  1  FUEL  LOADINGS,  5  AREAS,  AND  3  TYPES  OF  FIRE  SPREAD.  7 

C  THE  FOLLOWING  ARE  INPUT  VARIABLES:  8 

C  EF,  ENERGY  CONTENT  OF  FUEL  (BTU/POUND);  DTV ,  VIOLENT  BURNING  TIME  9 

C  (SECONDS);  EFV,  VIOLENT  BURNING  FRACTION?  DTR,  RESIDUAL  BURNING  TIME  10 

C  (SFCONDS);  EFR,  RESIDUAL  BURNING  FRACTION;  AK,  RATE  OF  AREA  CHANGE,  11 

C  (SQUARE  FEET/SEC);  AL,  RATE  OF  PERIMETER  CHANGE  (FEET/SEC) ?  AM,  12 

C  RATE  OF  RAOTUS  CHANGE  (FEET/SEC)?  FL(3),  FUEL  LOADINGS  (TONS/ACRE);  13 

C  A ( 5 ) ,  AREAS  (SQUARE  MILES);  ND(8),  NUMBFR  OF  INITIAL  FIRES  PER  SQUAR  14 

C  MILE  (NUMBER/MILES  SQUARED);  T ( 101  *  TIMES  AFTER  IGNITION  (SECONDS).  15 

1  FORMAT  (8F10.4)  16 

2  FORMAT  (818)  17 

3  FORMAT  (3F10.6)  18 

READ  (1,1)  EF,OTV,OTR,FFV,EFR  19 

READ  (1,3)  AK, AL,AM  20 

READ  (1,1)  (FL( I ), 1=1,3)  21 

READ  (1,1)  (AU), 1  =  1,5)  22 

READ  (1,2)  ( ND( I), 1=1,8)  23 

READ  (1,1)  (T( I ) , 1  =  1, 10)  24 

AK2=(AL*AL)/6. 28318  25 

AK3=( AM*AM)*6. 28318  26 

FMF=2.78784F+07  27 

FDV=EFV/DTV  28 

EDR=EFR/DTR  29 

OTVR=  DT V+DTR  30 

C  ICASE  INDICATES  WHICH  TYPE  OF  FIRE  SPREAD  IS  BEING  USED:  31 

C  1  =  CONSTANT  AREA  CHANGE;  2  =  CONSTANT  PERIMETER  CHANGE;  32 

C  3  =  CONSTANT  RADIAL  CHANGF.  33 

DO  45  ICASF=l,3  34 

IF  ( I  CASE-2  )  6,4,5  35 

4  CON$T=AK?  36 

GO  Tn  6  37 

5  C0NST=AK3  38 

C  M  INDEXES  OVER  THE  3  EUFL  LOADINGS.  L  INDEXES  OVER  THE  5  AREAS  39 

6  DC  45  M=  1,3  40 

no  45  L  =  I , 5  41 

on  7  1=1,8  4? 

N I  I ) =  N n I  I )*A(L )  43 

.N( I )=N( I )**< 1 .0/3.0)  44 

AA=(  (EL(M)*Ef  )♦♦(  1. 0/3.0)  )/(  AU)**(  l.  0/6.0)  14.031501  45 

STATEMENTS  8-15  COMPUTE  VELOCITIES  FOR  ICASE=l.  46 

IF  (ICASF-l)  8,8,16  47 

8  R=  AK  48 

DO  15  K= 1 ,10  49 

T I MF*T ( K l  50 

IF  (TIMF-DTV)  9,9,10  51 

9  ASV*R*T  I**E  52 
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ASR=0.0  53; 

ASB0=0.0  54 

cn  TO  13  55; 

10  IF  ( T IME-DT VR )  11,11,1?  56 

11  ASV=R*DTV  57 

ASR=R*l T IME-DT V )  58 

ASB0=0.0  59 

GO  TO  13  60 

12  ASV=R*DTV  61 

ASR=R*OTR  62 

ASBQ=R* ( T IME-OTVR )  63 

13  AAA=AA*( <ASV*EDV*ASR*EDR)**( 1.0/3. Q) )  64 

DO  14  1=1,8  65 

14  VW(K, I)=AAA*CN(I)  66 

15  NC ( K ) = A ( L ) / ( AK^T IMF ) *FMF  67 

GO  TO  ?4  68 

C  STATEMENTS  16-23  COMPUTE  VELOCITIES  FOR  ICASE=2  OR  3.  69 

16  DO  23  K=1 , 10  70 

TIME=T(K)  71 

R=CONST*TIME  72 

RH=R/2.0*  73 

IF  (TIMF-OTV)  17,17,18  74 

17  ASV=PH*TIMF  75 

ASR=0.0  76 

ASB0=0.0  77 

GO  TO  21  78 

18  IF  ( T IME-OTVR )  19,19,20  79 

19  ASV=RH*(DTV*(?.0-(DTV/TIME» ) )  80 

ASR=RH*(TIME-DTV*(2.0-0TV/TIME) )  81  5 

ASB0=0.0  82 

GO  TO  21  83 

20  ASV=RH*OTV*( 2.0-DTV/TIMEJ  84 

ASR=RH*0TR*(?.0-(2.0*DTV«-DTR  )/TIME)  85 

ASRO=RH* (T IME-DT VR*( 2.0-OTVR /TI ME ) )  86 

21  AAA=AAM ( ASV*FDV+ASR*EDR ) *♦( 1 .0/3.0) )  87 

00  2?  1*1,8  '  88 

22  VW  ( K  » I )=AAA*CN( I )  89 

23  NC(K)=(2.0*A(L ) )/(CONST*TIME*TIME )*FMF  90 

24  WRITE  13,25)  91 

25  FORMAT  Cl',T39,'VW  IS  THE  RADIAL  INDRAFT  WIND  VELOCITY  IN  MILFS  P  92 

*ER  HOUR •/'O' ,T23, *N  IS  THF  NUMBER  OF  INITIAL  FIRES,  NC  IS  THE  REQH  93 

1IRED  N  FOR  COALESCENCE  BY  THE  GIVEN  TIME'/)  94 

C=DTV/60.0  95 

B=0TR/60.0  96 

WRITE  (3,26)  FF ,EFV,EFR,C,B  97 

26  FORMAT  <• 0 ', T41 ,' ENERGY  CONTENT  OF  THE  FUFL  IS  ',F9.l.'  BTU  PFR  PO  98 

tUNP'//'0', T9, 'VIOLENT  BURNING  FRACTION  IS  *»F6.?,T85, ' RESIDUAL  R'J°  99 

INING  FRACTION  IS  • , F6. 2/ 1 X, T9, • VI OLENT  BURNING  TIME  IS  ',F6.2,'  MI  I  00  t 
SNUTFS' ,T8S, 'RESIDUAL  RURNING  T I MF  IS  *,F6.2,'  MINUTES'/)  1 0 1  | 

IF  ( I  CASE-2 )  27,28,29  10?  \ 

21  WRITE  (3, 30)  AK  10)  j 

GO  TO  33  104  ] 
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28 


WRITE  (3,31)  AL  105 

GO  TO  33  106 

29  WRITE  (3,32)  AW  107 

30  FORMAT  ( • O' , T32, 'CASE  1,  CONSTANT  AREA  INCREASE,  K  =  *,613.6,'  (SO  108 

SUARE  FEET/SECOND)'//)  109 

31  FORMAT  ( *  O' , T33, 'CASE  2A,  CONSTANT  PERIMETER  INCREASE,  L  =  '  110 

$,  E13.6, '  (FEET/SECOND)*//)  Hi 

32  FORMAT  ( • O' , T34, 'CASE  28,  CONSTANT  RADIAL  INCREASE,  M  =  *,E13.6,'  112 

i( FEET /SECOND) • // 1  113 

33  WRITE  (3,34)  A( L I , FL ( Ml  114 

34  FORMAT  ( • O' , T26, • AREA  =  *,F10.4,'  SQUARE  MILES  ',10X,»  FUEL  LOADIN  115 

$G  =  * , F10.4, '  TONS  PER  ACRE'//)  116 

DO  35  1*1,10  in 

35  TM( I )*T ( I ) /60.0  11« 

DO  40  11*1,2  HO 

IF  (II-ll  36,36,37  120 

36  I  A*  1  121 

18*5  122 

GO  TO  38  123 

37  I A*6  .  124 

1 8=10  125 

38  WRITE  (3,41)  ( TM ( I ) , I  =  I  A,  I B )  126 

WRITE  (3,42)  127 

00  39  K=1 , 8  128 

39  WRITE  (3,43)  ( N( K ) , VW( I , K ) , I  =  1  A , I B )  129 

40  WRITE  (3,44)  ( NC ( I  )  , I* I  A,  18 )  130 

•*  41  FORMAT  (• O' ,T4,5(4X, •  FOR  ',F6.1,'  MINUTES  ')/)  131 

42  FORMAT  ( ' 0 • , T4, 5 ( 6X , • N • , 1 IX , • VW* , 4X ) // )  132 

;  43  FORMAT  ( 1 X, T4, 5 ( 3X , 1 6, 2X, El  3 . 6) )  133 

:  44  FORMAT  I • O' , T4, 5 ( 7X , *NC  = • , 11 0,3X ) //// )  134 

45  CONTINUE  135 

READ  (1,2)  N2  136 

IF  (N2-2)  10^,46,103  137 

C  PART  TWO  CALCULATES  THE  PARTIAL  DERIVATIVES  WITH  RESPECT  TO  THE  138 

C  FOLLOWING  VARIARLFS:  N, A , FL , EF , EF V, EFR , R.  FACH  PARTIAL  IS  EVALUATED  139 

C  FOR  ALL  THE  INPUT  VALUES  OF  THE  VARIARLE  AND  AT  FACH  OF  THE  TIMES  140 

C  USED  IN  PART  ONE.  HI 

C  ADDITIONAL  INPUT  NFFPFC  FOR  PART  TWO  IS:  14? 

C  NDR,AR,FLR,FFP,FFVR,FFRR,AMP  (THFSF  ARE  MUST  PROBABLE  OR  REPRESENT  AT  143 

C  VALUES  FOR  THE  VARIABLES  FOR  WHICH  THE  PART  I AL  S  ARE  EVALUATED.)  144 

C  E  F? ( 4 )  ,  FNERGY  CONTENT  OF  THE  FUFL ( RTU/POUND) ;  FFV2I4I,  VIOLENT  145 

C  BURNING  FRACTION;  FFR2I4),  RESIDUAL  BURNING  FRACTION;  AM?(4),  PATE  146 

C  OF  RADIAL  INCRFASF(FFET/SFC  ).  147 

*  46  READ  (1,47)  NDR , AR , FLR , FF P , F F VP , FF R R , AMP  14H 

47  FORMAT  (|8,6F10.2)  149 

j'  READ  (1,48)  (EF2( I  1,1*1,41  150 

1  RF AD  (1,48)  (FFV2(  t), 1  =  1,4)  151 

RE AO  (1,48)  (FE«2(  I  )  ,1=1 ,4)  15? 

.  RE AO  (1,48)  (AM?( I  ), 1  =  1,4)  1*3 

|  48  FORMAT  (4F10.5)  1*4  5 

C  THROUGH  53  CALCULATES  QU  THE  PART  OF  THF  PAR  T I  AL  S  WHICH  !S  ONIY  16S  • 

C  T  IME-OFPFNOFNT  ■ 


I 
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AK2=AMR*AMR*6. 28318  157  j 

DO  53  K-  l , 10  158  ! 

T IM£=T ( K )  159 

RRH* AK2*T1 ME/2 • 0  160 

IF  IT IME-DTV)  49,49,50  161 

49  ASVR  ( K )-RRH*T IME  162 

ASRR (K 1=0.0  163 

GO  TO  53  164 

50  IF  ( T IMF-DTVR )  51,51,5?  165 

51  ASVR(K)=RRH*(DTV*(2.0-0TV/TTMF) )  166 

ASRR (K)=RRH*(TIME-DTV*( 2.0-DTV/TIME ) )  167 

GO  TO  53  168 

52  ASRR(K)=RRH*(0TR*(2.0-f <2.0*DTV+DTR)/TIMF) ))  169 

ASVR(K)=RRH*(0TV*(2.0-DTV/TIME)  )  170 

53  Q3 (K  )= ( ( ASVR ( K ) *EFVR /DTV ) ♦(  A  $R  RIK  ) *EFRR/DTR ))**(!• 0/3.0)  171 

C  PARTIALS  WITH  RESPECT  TO  N( NUMRFR  OF  INITIAL  FIRES)  FOR  8  N  172 

8B=.031501*( ( FLR*£FP )**( 1 ,0/3.0)/( 3.0*C AR**( 1 .0/6.0) >) i  173 

00  54  1=1,8  174 

BBB=BP/(KO(  I)*AR)**(2. 0/3.0)  175 

00  54  K=1 , 10  176 

54  PVW(  I,K)=8Bfl*03(K)  177 

55  FORMAT  I » 1 • , T40 , • PART  I AL  DERIVATIVES  OF  THE  RADIAL  INDRAFT  WIND  VF  178 

HOCITYM  179 

56  FORMAT  (IX, T 44, 'WITH  RESPECT  TO  N ,  THE  NUMBER  OF  INITIAL  FIRES')  180 

57  FORMAT  (1X,T31,'F0R  CASE  28,  CONSTANT  RADIAL  INCREASE,  M  =  '  181 

S, E13.6, '  FEET  PER  SECOND'//)  182 

58  FORMAT  ( • O' , T4, • WITH  AREA  *  ',F10.5,'  SQUARE  MILES  FUEL  LOADIN  183 

$G  =  ' , F10. 4, '  TONS  PER  ACRE  ENERGY  OF  FUEL  =  *,F10.4,'  BTU  PER  184 

$  POUND* )  185 

59  FORMAT  ( IX, T9, • VIOLENT  BURNING  FRACTION  =  • ,F6. 2 , T85 , • RES  I  DUAL  BUR  186 

$NING  FRACTION  =  *,F6.2)  187 

60  FORMAT  (  1 X, T9, • V IOLFNT  BURNING  TIME  =  ',F6.2,*  MI  NOTES ' ,T 85 , ' RES  1 0  188 

SUAL  BURNING  TIME  =  ',F6,?,»  MINUTES'////)  189 

61  FORMAT  C0',T15,'N  ='  ,  8(  3X,  1 8,3X) )  190 

62  FORMAT  ( ' 0* , T3, F6. I, '  MINUTES  ' ,8( IX, E13.6) )  191 

WRITE  (3,55)  19 2 

WRITF  (3,56)  193 

WRITE  (3,57)  AMR  194 

WRITF  (3,58)  AR,FLR,EFP  195 

WRITF  (3,59)  FFVR.FFRR  196 

WRITE  (3,60)  C,B  197 

00  63  1=  i , 8  198 

63  NA ( I ) =N0 ( I ) *AP  199 

WRITF  (3,61)  (NA( I), 1=1,8)  200 

DO  64  K= l » 10  >01 

64  WRITF  (3,6?)  TM ( K ) , ( PVW ( l , K ) , I = 1 , 8 )  20? 

C  PARTIALS  WITH  RESPECT  TO  A( AREA)  FOR  5  A  203 

RB=. 031501* (FLR*EFP)**( 1.0/ 3.0)  204 

DO  65  1*1,5  205 

H8B=BB* ( NOR*AR ) *• C !. 0/3.0)/ ( -6.0*A (1) ♦*(7.0/6. 01)  206 

no  65  K*l, 10  >07 

65  PVW( 1 ,K)*BBB*Q3IK)  >08  f 


f 
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NR  =  N[)R*AR  200 

WRITE  (3,55)  210 

WRITE  (3,66)  211 

66  FORMAT  ( IX, T51, • WITH  RESPECT  TO  A,  THE  TOTAL  AREA')  212 

WRITE  (3,57)  AMR  213 

WRITE  (3,67)  NR,FLR,FFP  214 

67  FORMAT  ( •  0 '  ,  T4,  •  W I  TH  INITIAL  NUMBER  OF  FIRFS  =  '  ,  I  8  ,.3X , 'FUEL  LOAOI  215 

$NG  =  * , F10.4,  '  TONS  PER  ACRE ' ,  3X» 'ENERGY  OF  FUEL  =  ',F10.4,'  BTU  P  216 
$FR  POUND' )  217 

WRITE  (3,59)  FFVR, 6FRR  218 

WRITF  (3,*0)  C,B  219 

WPITE  (3,69)  (All), 1=1, 5)  220 

on  68  K= 1 , 10  221 

68  WRITF  (3,70)  TM(K),(PVW(I »K) ,1=1,5)  222 

69  FORMAT  C0',T19,'A  *  •  ,  5  (  8X ,  F 10 . 5 , 2X)  /  IX  ,  Tl  3 ,  •  (SOUARE  MILES)'//)  223 

70  FORMAT  ( ' O' , T7, f 6. I, •  MINUTES  • , 5 ( 7X , E 13 . 6 ) )  224 

C  PARTI ALS  WITH  RESPECT  TO  FL ( FUEL  LOADING)  FOR  3  FL  225 

BB=.031501*(NDR*AR*EFP)**( 1 .0/3.0) /( 3.0*AR**( I. 0/6.0) )  226 

DO  71  1=1,3  227 

RRB=BB/FL(  I  )**( 2. 0/3.0)  228 

no  71  K  = i ,  10  229 

71  PVW( I , K )  =  BBB*03 (K  )  230 

WRITE  ( 3,56)  231 

WRITE  (3,72)  232 

72  FORMAT  ( IX , T49, » W I TH  RFSPFCT  TO  EL,  THE  FUEL  LOADING')  233 

WRITF  (3,57)  AMR  234 

WRITE  (3,73)  NR , AR , FFP  235 

73  FORMAT  ( • 0 • , T4 , • W I TH  INITIAL  NUMBER  OF  FIRFS  =  ' , 1 8 , 5X , • ARF A  =  •  236 

FfFlO.5, *  SOUARF  M ILFS •, 5XV • ENERGY  DP  FUEL  =  *,F10.4,'  RTU  PER  POUN  ?37 
*n*)  238 

WPTTF  (3,59)  EFVR.EFRP  239 

.  WRITF  (3,60)  C , B  240 

WRITF  (3,75)  (FL(  I  >,1  =  1,3)  241 

DU  74  K=l, 10  24? 

74  WRITF  (3,76)  TM ( K ) , ( PVW ( I ,K ) , I = 1 , 3 )  243 

75  FOPMAT  ( • 0 • , T 3? ,  ' FL  = • , 3 ( 1 1 X , F 1 0. 4 , 2X ) / 1 X , T26 , ' ( TONS  PFR  ACKF  ) »// )  244 

76  FORMAT  (*0',T2?,F6.l,'  M I NU TF S • , 3 ( l OX , F 1 3 . 6 ) )  245 

C  PARTIALS  WITH  RFSPFCT  Tn  FF ( ENERGY  CONTFNT  OF  FUEL)  FOR  4  EF  246 

9R=  I  •  f)  3  1  50 1  *  (  NOR*  A  R  *FL  R  )**(  l. 0/3.0)  5/(3. 0*AR**(  l  .0/6.0)  )  ?4  7 

DD  77  1=1,4  248 

B9B=RR/FF2( I )**(2. 0/3.0)  249 

DO  77  K=l,  10  250 

77  PVW( I,K)=RBR*03(K)  251 

012=.03 1501* (NOR *AR*FLP*FFP)**( l . 0/ 3. 0) / AR **( l . 0/6. 0 )  ?5  2 

WPITF  (3,55)  253 

WRITE  (3,78)  264 

78  FORMAT  (  IX,  T42, 'WITH  RFSPFCT  TO  FF  ,  THF  'ENERGY  CONTENT  OF  THF  FUEL  255 

J'>  254 

WRITE  ( 3,57)  AM°  257 

WRITF  (3,79)  NR , AR ,FLP  25H 

79  FORMAT  ( 'O' ,  T4,  •  WITH  INITIAL  NUMBER  OF  FIRES  =  •  ,  I  8 , 5X  ,  •  ARF  A  •  ?59 

ttE10.5, •  SQUARE  MILFS' ,5X, 'FUEL  LOAOI NG  =  *,F10.4,'  TONS  PF«  AC»F'  260 


III- 8 


i)  261 

WRITE  (3,59)  EFVR ,  EFRP  26? 

WRITE  (3,60)  C , R  263 

WRITE  (3,80)  (EF2( I ), 1=1,4)  264 

80  FORMAT  (»0',T?5,'FF  = ' , 4( 1 IX , F 10, 5 , 2X ) / 1 X ,T17 ,  •  ( BTU  PER  POUND)’//)  ?65 

DO  81  K  =  1 , 10  266 

81  WRITE  (3,82)  TM(K),(PVW(I,K), 1=1,4)  267 

82  FORMAT  ( • 0* , T14, F6. I, •  MINUTFS  • ,4( 10X,E13.6) )  268 

C  PARTIALS  WITH  RESPECT  TO  FFV( VIOLENT  BURNING  FRACTION)  FOR  4  EFV?  269 

DO  83  1=1,4  270 

EFRC=1.0-EFV?(I)  ?71 

DO  83  K= l , 10  272 

AAA= ASVR (K ) /DTV  273 

BBB=ASRR(K)/DTR  274 

8  3  PVW(I,K)  =  012*( AAA-BBB )/(3,0*( AAA*EFV2 (  I )  +  BBB*EFRC ) ** ( 2. 0/  3 .0 )  )  275 

WRITE  (3,55)  276 

WRITE  (3,84)  277 

84  FORMAT  ( IX, T42, *  WITH  RESPECT  TO  EFV,  THE  VIOLENT  BURNING  FRACTION*  278 

*)  *  279 

WRITF  (3,57)  AMR  280 

WRITE  (3,79)  NR,AR,FLR  281 

WRITE  (3,85)  EFP  28? 

85  FORMAT  ( IX, T9, ' FNERGY  OF  FUEL  =  *,F10.4,'  BTU  PER  POUND  * ,TR5, *RES I  283 

SDUAL  RURNING  FRACTION  =  1.0-EFV*)  284 

WRITE  (3,60)  C , B  295 

WRITF  (3,87)  ( FFV2( 11,1=1,4)  296 

00  86  K*l,10  287 

86  WRITE  (3,82)  TM( K )  ,  ( PVW( I ,K ) , 1  =  1 ,4)  288 

87  FORMAT  ( *0* ,T24, *EFV  = • , 4( 1 1 X , F 10.5.2X l // )  ?R9 

C  PARTIALS  WITH  RESPECT  TO  FFR(P6$IDUAL  BURNING  FRACTION)  FOR  4  EF9?  290 

DO  88  1=1,4  291 

EFVC=1.0-FFR?(I  )  292 

00  88  K= 1 , 10  293 

BRB=ASRR(K)/DTR  ?94 

A A A= ASVR ( K ) /DT V  295 

88  PVW( I,X)=012*(BBB-AAA)/(3,0*(AAA*EFVC+BBR*EFR2(l ))♦*(?. 0/ 3 .0 ) )  296 

WBITF  (3,55)  297 

WRITF  (3,89)  298 

89  FORMAT  ( 1X»T42, 'WITH  RESPECT  TO  EFR,  THE  RESIDUAL  BURNING  FRACTION  29** 

$')  300 

WRITF  (3,57)  AMR  301 

WBITF  (3,7**)  NR,AR,FLR  30? 

WRITE  (3,90)  FFP  303 

90  FORMAT  ( IX, T9, 'ENERGY  OF  FUEL  =*,F19.4,'  BTU  PER  POUND  * ,T85,*VIOLF  304 

TNT  BURNING  FRACTION  *  1.0-EFR')  308 

WRITE  (3,60)  C,R  306 

WRITF  ( 3, °2 )  (EFR2( I  ), 1  =  1,4)  307 

no  s»i  k*i, io 

91  WRITF  (3,82)  TM(K),(PVW(  I  ,K)  ,1*1,4)  MV7 

9?  FORMAT  (*0,,T24,'EFR  «',4(  U X,F 10.5, ?X ) / / )  3)0 

C  PARTIALS  WITH  RFSPFCT  TO  ®(CASE  ?B)  FOR  4  AM  Ml 

0123=QI?*EFVR/0TV  3)2 
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on  97  1*1,4  313 

CONST=AM2(I)*AM2( I )*6.28318  314 

00  97  Ks l ,  10  315 

TIMF*T(K)  314 

R=C0NST*TIMF/2.0  317 

IF  ( T IME-DTV )  93,93,94  318 

93  PVW(I,K)=Q123*TIM£/(6.0*(R*TIME*EFVR/DTV)**(2. 0/3.0))  319 

GO  TO  97  320 

94  IF  ( T IME-DTVR )  95,95,96  321 

95  8B*Q1 2* ( ( 2.0-OTV/T IME ) *EFVR*( (  T1ME-DTV*{ 2 .O-DTV/TIME) ) /OT  R*EFRR) )  322 

QQ*(R«(DTV*(2. O-DTV/TIME)  )*EFVR/DTV)*(  R*(  TIME-DTV*  I  2.0-DTV/T  IME >  »  323 

$*EFRR/DTR )  *24 

PVW(I,K)»BB/(6.0*QQ**(2.0/3.0) )  325 

GO  TO  97  326 

96  RB=Q12*( <( 2. O-DTV/TIME )*FFVR)  +  <  <2.0- (OTVR+DTV) /TIME) *FFRR)  )  327 

QQ= (R*DTV*(2.0-0TV/TIMF)*EFVR/DTV)+(R*(DTR*(2.O-( 2. 0*DT V+DTR )  328 

*/TIME) JOEFRR/OTR)  329 

PVW (I»K)*BB/( 6.0*0Q** (2.0/3. 0) )  330 

97  CONTINUE  331 

WRITE  (3., 55)  332 

WRITE  (3,98)  333 

98  FORMAT  ( 1X,T26, ’WITH  RESPECT  TO  R,  THF  RATE  OF  SPREAO,  FOR  CASE  2B  334 

$,R  *  6.28*M**2»TIMF*  )  335 

WRITE  (3,99)  AR ,NR ,FLR ,EFP  336 

99  FORMAT  (  •  0 • , T  14,  •  W I TH  AREA  *  ',F10.5,'  SQUARE  MILES', T75,' NUMB FR  0  337 

$F  INITIAL  FIRES  =  •  ,  I8/1X.T19, 'FUEL  LOADING  =  ',F10.4,*  TONS  PER  A  338 
$CRF' ,T75, ’ENERGY  OF  FUFL  =  *,F10.5,’  BTU  PER  POUND')  339 

WRITE  (3,100)  FFVR , EFRR ,C , B  340 

100  FORMAT  (  IX, T19, 'VIOLENT  BURNING  FRACTION  =  • ,F6. 2 , T75 , • RES  I  DU AL  RU  341 

$RNING  FRACTION  *  • .F6.2/IX, T19, 'VIOLENT  BURNING  TIME  =  ',F6.2,'  MI  342 

$NUTES  *,T75, 'RESIDUAL  RURNING  TIME  =  ',F6.2,»  MINUTES'////)  343 

WRITE  (3,101)  ( AM2 ( I), 1=1, 4)  344 

101  FORMAT  ( ' 0 • , T26, *M  *',4(10X,F13.6)/1X,T18,' (FEET /SECOND ) ' / /  )  345 

DO  102  K*l,10  346 

102  WRITE  (3,82)  TM ( K ) , ( PVW( I, K), 1=1, 4)  347 

103  CONTINUE  348 

STOP  349 

END  350 
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PARTIAL  DERIVATIVES  Lip  THE;  RACIAL  INDRAFT  WINC  VELOC  ITT 
WITH  RESPECT  TO  A,  THE  TOTAL  AREA 

FUR  CASE  28,  CONSTANT  RACIAL  INCREASE,  M  =  O.IOOOOOE-Ol  FEET  PER  SECOND 
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APPENDIX  IV 


RADIAL  INRUSH  WIND  VELOCITY  CURVES 

The  following  45  sets  of  curves  are  based  upon  the  rationale  and 
parameter  values  discussed  in  Appendix  IA  and  are  drawn  from  the  output 
of  Part  I  of  the  computer  program  listed  in  Appendix  III. 

Each  curve  on  a  page  represents  radial  inrush  wind  velocity  (v  )  as 

a  function  of  the  number  of  initial  ignitions  (n)  for  a  specified  time  after 

ignition.  Each  family  of  curves  shows  ten  different  times  with  a  specified 

combination  of  spread  case,  area,  and  fuel  loading. 

* 

The  curves  are  funded  above  by  the  coalescence  locus,  as  discussed 
in  Appendix  ID. 

For  all  of  the  curves  in  this  group,  the  fuel  energy  is  taken  as 
7000  Btu/lb,  the  violent  burning  time  as  6  minutes,  the  residual  burning 
time  as  70  minutes,  the  violent  burning  fraction  as  0.4,  and  the  residual 
burning  fraction  as  0.6  . 
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APPENDIX  VI 


ACCELERATED  SPREAD  STUDY 

VI.  1  General 

Given  n  initial  point  fires  distributed  uniformly  in  an  area  A  such 
that  the  initial  separation  distances  between  any  fire  and  its  nearest  neigh¬ 
bors  are  equal,  then  D,  x,  and  r^  are  defined  in  Fig.  VI- 1  for  any  two 
fires. 


Fig.  VI-1 

Assume  that  r  =  r  =  d(t) ,  that  both  fires  maintain  essentially  circular 
c  c 

1  2 

fronts  as  the  fires  spread,  and  the  radius  at  any  time  t  is  r  for  each  fire. 

c 

Let  M  represent  the  rate  of  radial  spread  in  ft/ min ,  so  that 

dr 

M  =  .  (VI.  1) 

dt 

M,  x.  and  r  are  to  be  determined  as  functions  of  time  for  four  ini- 

c 

tial  spread  modes,  with  accelerated  spread  beginning  at  some  value  of  r  . 


VI- 1 


* 

VI.  2  Accelerated  Spread  Relations 

Assume  that  for  r  >  [  VA/4n'  -  3F]  ,  i.  e. ,  when  t  >  t  ,  spread 

0  2. 

is  accelerated  according  to  the  hyperbolic  relation 

Mx  =  3F  •  M 

a 

where 

A  is  the  fire  area  (sq  ft), 

n  is  the  total  number  of  fires  in  area  A , 

F  is  the  flame  height  (assume  F  =  3  ft) , 

t  is  the  time  at  which  accelerated  spread  begins, 

2 

x  is  the  distance  of  one  fire  front  from  the 

point  of  merging, 

2x  is  the  minimum  separation  between  fires  (ft) , 

and 

M  is  the  value  of  M  at  ta ,  which  is  the  start  of 

di  u 

the  accelerated  spread,  i.  e. , 

r  at  t  =  \l A / 4n'  -  3F 
c  a 

NOTE:  D=  V A [n  ;  therefore,  since  x  =  ?D-r  , 

-  c 

x  at  t  5  3F  5  9  ft 
a 


VI.  3  Initial  Spread  Mode  One 

For  initial  spread  mode  one,  take,  for  0  <  r^  <  [  N  A  /  4n'  -  3F]  , 


Provided  by  Mr.  Craig  C.  Chandler,  Division  of  Fire  Heseaz  ch,  U.  S. 
Forest  Service,  Washington,  D.  C. 
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M  =  0  if  t  <  1/6760  min 
1  27 

M  =  logl0(6760t>  for 


6760  <  1  ^  ta  ) 


(VI.  3) 


Then,  since 


r  (t) 
c 


m  -  1.  27 

c  '  ~  2ff(ln  10) 


j  log10(6760t)dt  , 


6760 


t[ln(6760t)-l]  + 

— 


and 


(VI.  4) 


(VI.  5) 


x(t)  =  ]D-  r  (t)  =  *  VA/n  '  -  r  (t) 
c  c 


(VI.  6) 


where  r  (t)  is  the  value  ot  r  at  time  t, 
c  c 


VI.  4  Initial  Spread  Mode  Two 


For  initial  ~ oread  mode  two,  take,  for  0  <  rc  <  f  v  A  / 4n  -  3F]  , 


M  =  0  if  t  <  1 /0. 01  min 

M  =  log  (0.  01  t)  for 

^  7T  1  U 


—  <  '  i 


(VI.  7) 


As  before,  integration  of  Kq.  (VI.  7)  yields 


Provided  by  Mr.  Craig  C.  Chandler,  Division  of  Fire  Research,  U.  S. 
Forest  Service.  Washington,  D.  C. 
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(VI.  8) 


M  -  6.  67 

V  *  '  2ir(lnlO) 


t[ln(0.  91 1)  -  1  ]  + 


0.  91 


and,  also 


x(t)  =  \  D  -  r  (t)  =  \  N/A/n'  -  r  (t) 


VI.  5  Initial  Spread  Mode  Three 

Consider  now  the  case  of  constant  rate  of  radial  increase  for 


t  <  t  ,  with 
—  a 


In  this  case, 


and,  as  before 


M  =  C.  01  ft/sec  =  0.  6  ft /min 


r  (t)  =  0.  6 1 
c 


x(t)  =  1 D  -  r  <t) 
c 


VI.  6  Initial  Spread  Mode  Four 

Also,  consider  the  case  of  constant  rate  of  area  increase  up  to 

t  =  t  ,  i.  e.  , 
a 


dA  •>  '> 

— —  =  K  =  5  ft" /sec  =  300  ft" /min 
dt  1 


Now 


a 

dUr")  dr 

dA  c  „  c  _  „ 

__  a  — —  =  >nr  .  2 nr  M 

dt  dt  c  dt  c 


and  equating  (VI.  13)  with  (VI.  14)  yields 


2jrr  M  =  300  ft" /min 
c 


(VI.  9) 

(VI.  10) 

(VI.  11) 

(VI.  12) 

(VI.  13) 

(VI.  14) 

(VI.  15) 
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Taking  Eq.  (VI.  15)  and  remembering  that 


dr 

C  1 

M  “  — =  r 
dt  c 


we  have  that 


r  r 
c  c 


300 
2  7T 


=  C 


with  the  initial  condition  that  r  (0)  =  0  . 

c 

Hence, 


r  (t)  =  ^/2C^  =  ' 

C  V  7T 


and,  as  before 


x(t)  =  |  D  -  r  (t) 

c 


Also,  differentiation  of  Eq.  (VI.  17)  yields 


M  = 


75 


7T  t 


(VI.  16) 


(VI.  1 7) 


(VI.  18) 


(VI.  19) 


VI.  7  Method  of  Calculation 

A  computer  program  has  been  written  based  upon  the  foregoing 

which  calculates  M,  r  ,  and  x  as  functions  of  time  for  0  <  r  <  ’  D.  A 

c  —  c 

listing  of  this  program  and  its  input  requirements  is  included  in  Appendix  VII, 

VI.  8  Results 

The  results  of  calculations  with  the  above  program  are  presented 
graphically  in  Appendix  VIII. 
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APPENDIX  VII 


ACCELERATED  SPREAD  PROGRAM 


VII.  1  General 

The  following  program  is  based  upon  the  accelerated  spread  rationale 
explained  in  Appendix  VI,  and  calculates  spread  rate  (M) ,  fire  radius  ( r  ) , 
and  distance  to  merging  (x)  for  each  of  the  four  initial  spread  modes  dis¬ 
cussed  in  Appendix  VI,  with  acceleration  for  each  mode  beginning  when 
x  <  3F ,  where  F  represents  the  flame  height.  The  program  is  written  in 
FORTRAN  IV  for  an  IBM  System/  360  Model  40G  digital  computer.  The  pro¬ 
gram  is  in  three  parts.  The  first  part  performs  the  calculations  for  each 
initial  spread  mode.  SUBROUTINE  PRINT  prints  the  results  of  these  cal¬ 
culations  up  to  the  point  of  acceleration,  and  SUBROUTINE  ACC  performs 
the  calculations  and  prints  the  results  for  each  initial  mode  after  accelera¬ 
tion  begins. 

VII.  2  Data  Cards 

The  first  card  has  format  112  and  contains  NSETS,  the  number  of 
data  sets  to  be  used  in  the  run. 

Each  data  set  consists  of  a  single  card  with  format  5E12.  4  contain¬ 
ing  the  following  inputs: 


VII- 1 


Column 

Program 

Symbol 

Units 

Variable 

1-12 

D 

feet 

Initial  separation  between  fires 

13-24 

F 

feet 

Flame  height 

25-36 

TINC 

minutes 

Time  increment  of  calculation 
for  initial  spread 

37-48 

TM1 

minutes 

Spread  start  time  for  Mode  1 

49-60 

TM2 

minutes 

Spread  start  time  for  Mode  2 

NOTE:  Due  to  the  form  of  the  equations  for  initial  spread 
modes  one  and  two,  TM1  must  be  greater  than  1/6760  min¬ 
utes  and  TM2  must  be  greater  than  1/0.  91  minutes. 


VII.  3  Program  Listing 

(See  following  pages.  ) 
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PROGRAM  SPREAD 

DIMFNSIUN  T(1500),M<  1500 ),RC( 1500),  X(  1500) 

REAL  M,MA 
P I  =  3 .  14159 
READ  (1,1)  NSETS 
FORMAT  (II?) 

NSFTS  IS  THE  NUMBER  UF  OATA  SFTS 

00  21  K=1 , NSET  S 

READ  (1,2)  D* F «  T INC • TM1 , TM2 

D  IS  THE  INITIAL  SEPARATION  BETWEEN  FIRES  IN  FEET 
F  IS  THE  FLAMF  HEIGHT  IN  FEFT 

TINC  IS  THE  TIME  INCREMENT  IN  MINUTES  FUR  INITIAL  SPREAO 

TM1  IS  THE  SPREAD  START  TIME  FUR  MODE  ONE 

TM2  IS  THF  SPRFAD  START  TIME  FOR  MODE  TWO 

RC  IS  THE  FIRE  RADIUS  IN  FFFT 

X  IS  THE  DISTANCE  TU  MERGIMG  IN  FEET 

M  IS  THE  RATE  OF  RADIAL  SPREAD  IN  FEET  PER  MINUTE 

T  IS  TIME  IN  MINUTES 

FURMAT  (5E12.4) 

MODE-O 
T (  1 )  =  T I NC 

THIS  LOOP  CALCULATFS  SPREAO  FOR  MODE  ONE 

TFST=3.0*F 

DO  7  1=1,1000 

IF  (T(I)-TMl)  3,3,4 

M ( I  )  =0 

RC ( I ) =0 

GO  TO  5 

M(I)  =  I. 2  / /(2. 0*PI)  *AL  Of*  10(6760. 0*T(I  )  ) 

PC ( I  )  =  l  •  27/  (  2  *  0*P  I  *AL()G( 10.0) )*(T( l ) * ( ALUG ( 6 7faO. 0*T ( I > 1-1.0) 
4*1. 0/6760.0) 

X(  I  )  =0 . 5*U-kC ( l ) 

I  1  =  1 

IF  (  X (  I  ) -T E ST  )  8 , H ,  6 

T(H-l)»(Ul)*TINC 

CONTINUF 

MA  =  M(  I  1 ) 

T  A  =  T (II) 

X  A=  X (  I  I  ) 

R A=RC (  I  I  ) 

M,,OE  =  MODE  4-1 

CALL  PRINT  (X,M,T,RC»II, MODE  ,0 , F ) 

CALL  ACC  (XA,MA,TA,RA,X,M,T,RC,Mt|.)F,F) 

(ill  TO  (9,16,18,21)  , MODE 
T( I l=T INC 

THIS  LOOP  CALCULATES  SPREAD  FOR  ML"*E  T *0 
DO  14  1=1,1000 
IF  (TID-TM?)  10,10,11 

M (I) a  n . 0 
RC ( I  1  =  0.) 
go  rr  i? 

M(  I  1  =  6.67/ (  2.n*Pl  )*ALU",  1 0  (  0.  5l*T(  1)1 
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l. 


12 


13 

14 

15 
C 


16 
1  7 
18 
C 


19 

?0 


5 

r 


^C(  1  1  =  6.67/ ULGul  10.0)*2.0«'PI)*(Tm*(ALOG(0.9l*TU  11-1.0) 

$+1.0/0.91 ) 

X(  I  )  =  0 . 5*U-RC  (  1  1 
11  =  1 

IF  (X(  I  1-TEST)  dfdf  13 
T( I+l»=( 1+1 )*T1NC 
CONTINUF 
T  ( 1  )*T INC 

THIS  LOUP  IS  FOR  SPREAD  MODE  THREE,  CONSTANT  RADIAL  INCREASE 

DO  17  1=1,1000 

M (  11=0.6 

RC ( I >  =  0.6*T ( I  1 

XI  I  )=0.5*D-RCl I  1 

I  1  =  1 

IF  ( X  C I  1-TFST 1  8,8,16 
T(I+11=(I+11*TINC 
CONTINUF 
T  (  1 1 =T INC 

THIS  LOOP  CALCULATES  SPREAD  FOR  MODE  FOUR 

DO  20  1=1,1000 

M ( I  )=SURT(75.0/(PI*T( I  1 1 1 

RC ( I >  =  SURT( (  300 . 0*  T ( I  1  1/PI  1 

X( I  1=0. b+O-RC ( I  1 

I  1=1 

IF  (XIIl-TtST)  rt,«,19 

Ttl+1)=( I+1)*TTNC 

CONTINUF 

CONTINUE 

RETURN 

END 

SUBROUTINE  °RINT  l  X ,  M,  T  ,  RC. ,  I  I  , MODE  ,  0  ,  F  1 
DIMENSION  X(l),M(ll,Tm,RC(l) 

FORMAT  (  1H1, 10X, 19HINITIAL  SFPARA T I QN= , F 1 2 . 4 , 10X , 1 3HFL AME  HEIGHT; 
$,F 12.4, 10X, 21HSPREAU  MODE  IS  NUMBER, 151 
NPG=l  1/96 
L=  I  I-NPG*96 
M 1  =  1 
M2  =  9fc 

00  4  K=  1 , NPG 

WRITF  (3,11  U  »  E  ,  MODE. 

W*ITE  (3,2) 

Format  ( l«u, ix,?( l2X,4HT!Mt ,5X,11HSPRFAD  RATE , 5X , 1 1 HF I RF  RAOIUS 
$, 3X, 14HMLPUL  DISTANCE!! 

WRITE  (3,31  (T ( J  1 , Ml J 1 ,RC(JJ,X(J1  ,J  =  Ml,M2) 

EORMAT  ( 2X , 8F 1 6 .4 1 

M  t  =  M  1  +  46 

**2=M  2  +  9b 

CONTINUE 

IF  (L 1  9,9, S 

IF  (NPG)  6,6,7 

91=1 

M2  =  L 


53i 

54  j 

55 

56  | 

57 

58  I 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 
87 
08 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102  ? 

103 

104 
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nn  tp  h  105 

7  M2=M1+L-1  106 

8  WRITE  (3,1)  D ,  F ,  MODE  107 

WRITE  (3,7)  108 

WRITE  (3,3)  (T( J),M( J ) , RC < J) ,X( J) ,J=M1,M2)  109 

n  RETURN  110 

FND  111 

SUBROUT 1NL  ACC  ( XA , MA , 1  A , BA , X , M, T , RC , MUDE , F )  112 

C  THIS  SUBROUTINE  CALCULATES  SPREAD  RATE,  RAOIUS,  AND  113 

C  SEPARATION  DISTANCE  AFTER  ACCELERATION  BEGINS  114 

DIMENSION  XC  l),T(l),M( 1),RC(1)  115 

REAL  M, MA  116 

DT- .005  117 

X ( 1 )  =  XA  118 

M ( l ) *  MA  119 

T ( 1 ) =T  A  120 

RC ( 1 )-RA  121 

DO  l  J  =  2 »  l 500  122 

X( J )  =  X  C J-l)-M( J-1)*DT  123 

M(J)=3.0*F*MA/X( J)  124 

RC ( J ) =R A+XA-X ( J )  125 

T ( J ) =T  A  +  DT* ( J- 1 )  126 

IF  (XU)-. 05)  2,2,1  127 

1  CONTINUE  128 

2  LPINC=4  129 

N= J-L  P I NC  + 1  130 

LPPG=96*LPINC  131 

NPG=N/LPPG  132 

L=N-NPG*LfPG  133 

*1  =  1  -  134 

M  2  =  L  P  P  G  135 

DO  3  KK= l , NPG  136 

WRITE  (3,8)  MCOE  137 

WR)TC  (  j,S)  13* 

WRITE  (  3,10)  ( T (K) ,M(K ) , RC ( K ) ,X(K) ,K=M1 ,M2,LPINC)  139 

M 1  =  M  1  +L  PPG  140 

M2  =  M?  *L  PPG  141 

3  CONTINUE  142 

IF  (L)  11,11,4  143 

4  IF  (NPG)  5,5,6  1^4 

5  M 1 = 1  145 

M2  =  L  146 

r,i:  ru  7  147 

6  M7-Ml*L-l  146 

7  WHITE  (3,8)  '•UOt  149 

WRITE  ( 3,5)  160 

WRITE  (3,10)  ( T(K) ,m(k ) ,RC(K) ,X(K) ,K=Ml ,M2,LPINC)  151 

H  c  t.  P  M  A 1  (  l),l  ,40X, 'ACCLLLRAT  ItIN  PHASE  FOR  INITIAL  SPREAD  MODE', 15)  152 

9  FORMAT  ( 1H0, IX, 2( IPX, 4HTIMF ,5X, l IHSPREAD  RATE ,5X,1 IHFIRE  RAOIUS  153 

i, 3X, I4HFFf GE  DISTANCE))  154 

10  f  OK  MAT  (7X,5F16.4|  155 

11  RETURN  166 
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APPENDIX  VHI 


ACCELERATED  SPREAD  CURVES 


The  following  eight  curves  show  spread  rate  (M)  as  a  function  of 
time  (t)  and  fire  radius  (rc)  for  each  of  the  four  initial  spread  modes  dis¬ 
cussed  in  Appendix  VI,  with  the  point  at  which  accelerated  spread  begins 
indicated  by  the  symbol  A  beside  each  curve. 

These  curves  are  drawn  from  the  output  of  the  computer  program 
described  in  Appendix  VII,  with  an  initial  separation  between  fires  (D)  of 
78.  8  feet,  a  flame  height  (F)  of  3.  0  feet,  TINC  =  0.  1  minutes, 

TM1  =  .  000148  minutes,  and  TM2  =  1.1  minutes. 

The  value  of  D  was  taken  from  the  formula  D  =  \l A/n'  ,  with 
n  =  A*  .  Using  the  representative  values  of  A  and  n^  as  discussed  in 
Appendix  II , 


n/W  =  78.  8  ft 
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